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摘 要 : 基于 1998—2013 年 的 SEOTvec NDVI 数 据 ,利用 像 元 二 分 法 务 


板 校正 标准 误 (PCSE ) 估计 方法 ,分 析 陕 


西 省 植被 覆盖 时 空 变化 特征 以 及 退耕 还 林 工 程 .降雨 水 平等 自然 .社会 经 济 因素 对 植被 覆盖 变化 的 相对 贡献 ,探究 
退耕 还 林 工 程 对 植被 恢复 的 驱动 效应 。 结 果 表 明 :GD 陕西 省 植被 覆盖 总 体 呈 增加 趋势 ,1998 一 2013 年 年 均 归 一 化 


植被 指数 (NDVI) 值 增长 0.98% , 


在 退耕 还 林 工 程 实施 阶段 (1998 一 2009 年 ) ,植被 恢复 效果 更 为 明显 。@) 陕西 


省 植被 恢复 效果 空间 分 异 显著 ,16 a 间 , 各 地 区 NDVI 值 上 升幅 度 不 同 ,40. 61% 的 区 域 提 升 约 0.1,29.05% 的 区 域 
提升 约 0. 2,23. 69% 的 区 域 提升 0.25 以 上 。 其 中 ,位 于 陕 北 退 耕 还 林 区 的 延安 市 植被 恢复 效果 最 为 明显 ,NDVI 年 


增长 约 为 1.12% 。@ 退耕 还 林 工 程 对 植被 恢复 影响 显著 ,退耕 还 林 面 积 占 土 地 总 面积 每 提升 1% 将 使 NDVI 值 增 
加 1.97% ,相当 于 降水 量 增加 了 4.29% ,这 对 于 降水 量 少 且 不 均 的 陕西 省 生态 环境 建设 具有 重要 参考 意义 。 
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当前 ,中 国生 态 文明 建设 正 处 于 快速 推进 的 重 
要 时 期 ,如何 有 效 解决 由 以 往 不 尽 合理 的 经 济 发 展 
模式 而 导致 的 土地 沙化 .水 土 流失 等 生态 环境 问题 
已 经 成 为 学 术 界 及 各 级 政府 部 门 关注 的 重点 。 植 被 
作为 生态 系统 的 主体 部 分 " ,是 指示 生态 环境 变 
迁 的 重要 指标 4-53 。 植 被 恢复 是 生态 系统 恢复 的 
首要 工作 名。 退耕 还 林 工 程 作为 我 国 乃 至 世界 上 
规模 最 大 的 生态 修复 和 生态 补偿 工程 ,是 应 对 环境 
变化 和 人 类 活动 干扰 的 重要 措施 之 一 。 随 着 工程 的 
深入 推进 ,西部 退耕 区 植被 覆盖 度 显著 提高 "9 。 
如 何 从 区 域 尺 度 定量 分 析 植 被 恢复 进度 ,揭示 退耕 
还 林 工 程 对 植被 恢复 的 相对 贡献 ,对 指导 工程 顺利 
实施 以 及 为 深化 生态 文明 体制 改革 提供 决策 参考 具 
有 重要 现实 意义 。 


VD) 是 目前 广泛 应 用 于 反映 植被 状况 的 指标 , 它 与 
植被 覆盖 度 、 生 长 状况 .生物 量 和 叶 面 积 指 数 密切 相 
KOD ,能 够 在 较 大 时 空 尺度 上 反映 区 域 植被 覆盖 
信息 所] 。 学 者 们 从 国家 及 区 域 尺度 出 发 呈 -5 , 利 
用 NDVI 对 植被 覆盖 时 空 变化 特征 及 其 与 人 类 活动 
和 气候 变化 的 相关 关系 等 开展 了 大 量 研究 … O ,并 
取得 了 较为 丰富 的 成 果 和 可 资 借鉴 的 经 验 。 与 此 同 
时 ,部 分 学 者 重点 关注 了 退耕 区 的 植被 覆盖 变化 情 
况 ,信息 保 等 2 指出 , 自 退耕 还 林政 策 大 规模 实施 
以 来 ,黄土 高 原 植被 覆盖 度 显 著 提 高 ; 张 清 雨 等 2 
分 析 了 内 蒙古 不 同 生 态 区 内 NDVI 时 空 变 化 ,指出 
人 类 活动 对 NDVI 的 影响 程度 逐渐 增强 。 类 似 研 究 
很 多 ,结果 均 表 明 植 被 覆盖 的 时 空 变化 与 退耕 还 林 
工程 的 实施 密切 相关 。 然 而 ,上 述 研究 多 从 定性 角 


植被 覆盖 及 其 变化 是 区 域 生 态 环境 变化 的 重要 
指标 ” ,是 自然 人 文 等 因素 共同 作用 的 结果 ,尤其 
是 在 人 类 活动 影响 日 趋 强 烈 的 今天 ,植被 覆盖 变化 
更 能 深切 反映 人 类 活动 的 进程 。 近 年 来 , 随 着 地 理 
言 息 科学 的 发 展 与 成 熟 ,遥感 技术 已 成 为 获取 植被 
覆盖 相关 数据 的 主流 方式 。 归 一 化 植被 指数 ( ND- 
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度 描述 退耕 区 植被 变化 情况 ,从 定量 角度 评估 退耕 
还 林 工 程 对 植被 恢复 相对 贡献 的 研究 较 少 ,不 利于 
准确 评估 退耕 还 林 工 程 所 取得 的 实效 。 

陕西 省 是 我 国 水 土 流失 最 为 严重 的 省 份 之 一 ， 
省 内 地 貌 状况 复杂 ,南北 空间 差异 大 ,植被 覆盖 度 
低 ,整体 生态 环境 相对 脆弱 。 陕 北 地 区 在 20 t 
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纪 末 和 21 世纪 初 频 繁 发 生 极端 恶劣 天 气 ,生态 失调 
等 环境 问题 尤为 突出 。 鉴 于 此 ,为 有 效 治理 和 解决 
严重 的 生态 环境 恶化 和 水 土 流失 问题 ,陕西 省 于 
1999 年 率先 实施 退耕 还 林 工 程 ,历经 十 几 年 ,植被 
履 盖 度 显著 提高 , 林 草 地 面积 占 比 不 断 增加 ,生态 环 
境 明显 好 转 。 因 此 ,有 效 揭示 退耕 还 林 工 程 对 植被 
恢复 的 驱动 作用 ,对 于 推进 陕西 省 生态 文明 建设 具 
有 重要 的 意义 。 

本 文 基于 RS 和 GIS 技术 ,利用 SPOT), NDVI 
数据 ,提取 陕西 省 1998 一 2013 年 地 表 植 被 覆盖 变化 
的 时 空 信息 ,分 析 研 究 区 植被 覆盖 的 年 际 变化 和 区 
域 差 异 特征 ,并 通过 PCSE 估计 方法 定量 分 析 退 耕 
还 林 工 程 .降雨 水 平 人 地 关系 ,经济 发 展 水 平 .农耕 
强度 等 自然 ,社会 经 济 因 素 对 植被 覆盖 变化 的 相对 
贡献 ,进而 探讨 退耕 还 林 工 程 对 植被 恢复 效果 的 驱 
动 效应 ,以 期 为 客观 评价 陕西 省 退耕 还 林 工 程 的 生 
态 效益 提供 科学 依据 ,为 西部 地 区 深入 开展 生态 文 
明 建 设 提供 决策 参考 。 


1 研究 区 概况 


陕西 省 位 于 我 国 西北 内 陆 (31°42’ ~ 39°35'N, 
105°29' ~ 111°15'E) ,平均 海拔 1 127 m, HJE TŽ 
长 状 , 按 照 地 貌 可 划分 为 黄土 高 原 区 、 秦 岭 山地 区 、 
关中 平原 区 等 ,其 中 ,黄土 高 原 区 占 全 省 面积 最 大 ， 
占 比 达 40% 以 上 。 全 省 地 跨 北 亚热带 、 暧 温 带 和 中 
温带 3 个 气候 带 ,南北 气候 差异 显著 ,气温 由 南 向 北 
逐渐 降低 ,年 平均 气温 在 7 ~16 CC 。 年 降水 量 北 少 
PAZ ,年 际 降水 量 差异 较 大 , 且 降 水 各 季 分 配 不 均 
衡 , 多 集中 在 每 年 的 7 一 8 月 。 陕 西 省 水 土 流 失 严 
重 , 是 国家 退耕 还 林 工程 重点 建设 省 份 之 一 。1999 
年 起 ,陕西 省 启动 并 实施 了 退耕 还 林 工 程 。 根 据 陕 
西 省 林业 勘查 设计 院 提 供 的 退耕 还 林 区 图 (图 1)， 
陕西 省 退耕 还 林 区 包括 陕 北 退耕 还 林 区 和 陕 南 退耕 
还 林 区 ,分 别 位 于 陕 北 黄土 高 原 丘 陵 的 沟 密 区 、 渭 北 
拓 区 和 陕 南 低 山 丘 陵 区 。 陕 西 省 退耕 还 林 区 以 草地 
和 耕地 为 主 ,分 别 占 退耕 还 林 总 面积 的 45. 9% 和 
44.1% ;林地 面积 占 退耕 还 林 总 面积 的 9.0% ;水 
Jak ,城乡 .工矿 .居民 用 地 和 未 利用 地 所 占 比 例 不 及 
1% 5 。 自 退耕 还 林 工 程 实施 以 来 ,生态 建设 取得 
了 明显 成 效 ,十 多 年 来 累计 完成 退耕 还 林 计 划 任 务 
2.44 x 10° hm ,其 中 退耕 地 还 林 1.02 x 10° hm’, 
家 累计 投入 补助 资金 2.79 x10" 元 。 
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Fig.1 Schematic map of the areas where the Crain for Creen 


Project is implemented in Shaanxi Province 
2 数据 与 研究 方法 


2.1 数据 来 源 与 数据 预 处 理 

收集 了 1998 一 2013 年 陕西 省 10 个 主要 行政 区 
(市 /区 )16 a 的 数据 ,所 构成 的 面板 数据 总 计 有 960 
个 观测 值 。 
2.1.1 数据 ”数据 来 源 于 SPOTVysc NDVI 数据 产 
tu (Flemish Institute for Technological Research, Vi- 
to, Belgium) ,该 数据 从 1998 年 4 月 开始 提供 ,时 
间 精 度 10 qd, 空间 精度 1 km。 回 顾 以 往 国内 及 陕 
西 省 的 相关 研究 ,8 月 是 陕西 省 植物 生长 最 茂盛 的 
时 间 点 ,NDVI 值 处 于 全 年 最 高 水 平 "。 因 此 , 研 
究 使 用 ENVI 5. 1 软件 进行 格式 转换 ,提取 出 
1998 一 2013 年 每 年 8 月 底 的 波段 影像 图 。 由 于 原 
台 影像 已 经 过 辐射 校正 和 大 气 校正 等 预 处 理 , 故 
可 直接 利用 Band Math 工具 对 遥感 影像 图 做 格式 
转换 ,将 其 原始 数据 格式 (0 ~ 255 ) 转换 为 研究 需 
要 的 NDVI 值 ( 取 值 范围 ( -1,1] )2 ,计算 公式 
为 :NDVI= DN x0. 004 -0.1, 其 中 , NDVI (AATF 
零 表 示 和 覆盖 地 物 是 绿色 植被 , 且 值 越 大 表明 绿色 
植被 数量 越 多 ,和 覆盖 度 越 高 ; 值 越 接近 于 零 表 明 绿 


成 


=> 


色 植 被 覆盖 度 越 低 ; 若 NDVI 值 为 负 值 则 表示 地 表 
履 盖 物 是 非 植被 地 物 , 例 如 裸露 地 表 ( 兰 石 、 沙 
土 ) 水域 . 雪 、 冰 或 者 云 履 盖 等 。 

2.1.2 陕西 省 10 个 地 级 市 的 统计 数据 ”降水 量 数 
据 来 源 于 国家 气象 科学 数据 共享 服务 平台 ; GDP, 
CPI 数据 来 源 于 《陕西 统计 年 鉴 》; 人口 数 收集 于 国 
家 统计 局 ;耕地 面积 .土地 面积 收集 于 陕西 林业 局 ; 
退耕 还 林 面 积 收集 于 国家 林业 局 退耕 还 林 ( 草 ) 工 
程 管理 办 公 室 数据 共享 平台 。 

2.2 研究 方法 

2.2.1 植被 覆盖 度 提取 方法 ”植被 指数 (VI) 是 从 
遥感 影像 中 抽象 出 来 的 对 地 表 植 被 活动 简单 .有效 
和 经 验 的 量度 , 它 能 够 反映 地 表 植 物 生 长 状况 ,覆盖 
情况 .生物 量 、 叶 面积 指数 、 光 合作 用 强度 等 生物 物 
理 特征 。 大 量 人 研究 表明 ,植被 指数 与 植被 覆盖 度 之 
间 存 在 显著 的 相关 关系 ,用 它 来 计算 植被 覆盖 度 是 
合适 的 下。 其 中 , 归 一 化 植被 指数 (NDVI) 对 植 


册 昆 等 :陕西 省 退耕 还 林 工 程 对 植被 恢复 的 效应 


基于 PCSE 修正 的 面板 数据 模型 


被 类 型 的 NDVI 和 NDVI a RW EMIT f°" ， 
实际 计算 过 程 中 ,分 别 用 生长 季 内 NDVI 最 大 值 和 
最 小 值 代替 NDVI,,, 和 NDVI , 则 植被 覆盖 度 公式 
如 下 : 


f. = (NDVI - NDVL,,)/A( NDVI,,,, - NDVI.) (3) 


该 方法 已 在 其 他 地 区 广泛 应 用 ,并 且 精 度 较 
高 。 如 陈 效 还 等 2 用 此 方法 估算 内 蒙古 地 区 典型 
草原 植被 覆盖 度 ,并 对 估算 值 与 观测 值 进行 拟 合 ， 
拟 合 系数 接近 于 1, 表明 估算 结果 可 信 ; 李 苗 苗 
等 20 采用 此 方法 对 密云 水 库 上 流 进 行 了 植被 覆盖 
度 提取 ,并 与 实测 值 进 行 相关 分 析 , 相关 系数 达 
99% ,认为 估算 结果 可 靠 性 较 高 ,可 适用 于 大 面积 
估算 。 

2.2.2 面板 数据 模型 与 方法 ”研究 拟 选 用 面板 数 
据 回 归 模 型 来 直观 展现 退耕 还 林 工 程 对 植被 恢复 的 
影响 。 使 用 面板 数据 ,可 以 克服 时 间 序 列 分 析 受 多 


被 的 生物 物理 特征 十 分 敏感 , 且 在 时 效 、 尺 度 方面 都 
有 具有 显著 优势 ,通常 被 用 来 进行 区 域 尺 度 的 植被 分 
类 和 植被 覆盖 度 研究 。 
研究 采用 像 元 二 分 模型 建立 NDVI 与 植被 覆盖 
度 二 者 之 间 的 转换 关系 ,提取 植被 覆盖 度 信息 。 
像 元 二 分 模型 基于 混合 像 元 分 解 理论 ,将 影像 数据 
中 的 绿色 植被 和 裸 土 信 息 进 行 线性 分 解 ,进而 分 离 
得 到 像 元 中 的 植被 丰 度 信息 , 反 演 植被 覆盖 度 。 该 
方法 能 够 在 较 大 程度 上 消除 土壤 等 背景 因素 的 影 
响 , 且 对 影像 辐射 校正 影响 不 敏感 ,计算 简便 ,已 被 
广泛 应 用 于 计算 植被 覆盖 度 ,并 取得 了 较为 理想 的 
结果 。 
模型 假设 每 个 像 元 的 NDVI 值 由 绿色 植被 所 贡 
献 的 信息 与 裸 土 部 分 所 贡献 的 信息 合成 ,NDYVI 值 
NDVI =NDVI,,,f, + NDVI (1 —f.) 


(1) 


式 中 :NDVI 为 纯 植被 像 元 的 NDVI 值 ;NDVI 为 
纯 裸 地 像 元 的 NDVI 值 ; 人 为 植被 覆盖 度 。 根 据 公 
式 (1) ,的 计算 公式 如 下 : 

f. = (NDVI - NDVI,,;,)/(NDVI,,. - NDVI) (2) 


理论 上 ,NDVI,,.44 F 1,NDVI #3 F 0. 但 受 环 
Tae Al Ze AY , NDVI, FA NDVI 存在 一 定 的 变 
异 。 目 前 不 同 研究 对 NDVL.. 和 NDVI 的 取 值 方法 
有 较 大 差异 ,研究 结合 多 数学 者 的 取 值 经 验 , 将 各 植 


重 共 线 性 的 困扰 ,有 利于 观察 个 体 差异 ,增加 多 变 
性 ,减少 共 线 性 ,提高 自由 度 和 有 效 性 ,减少 偏差 。 
回归 模型 可 以 清晰 地 观察 因 变量 关于 另 一 个 ( 些 ) 
自 变 量 的 具体 依赖 关系 。 

为 了 避免 遗漏 重要 变量 造成 虚假 回归 ,在 把 退 
耕 还 林 工 程 作为 影响 植被 恢复 的 主要 自 变 量 外 ,还 
选取 降雨 水 平 ` 人 地 关系 、 经 济 发 展 水 平 农耕 强度 
作为 自 变量 对 植被 恢复 进行 整体 分 析 。 将 植被 覆盖 
水 平 作为 因 变 量 ,并 设立 随机 扰动 项 jw。 为 确保 人 
计 结 果 的 有 效 性 ,对 回归 方程 进行 对 数 化 处 理 , 并 利 
用 Stata 12.0 对 面板 模型 进行 Hausman 检验 ,得 P 
值 小 于 0.05 显著 性 水 平 , 故 拒绝 建立 随机 效应 模型 
原 假 设 ,选择 建立 个 体 固定 效应 模型 。 模 型 设置 如 
pO”, 

InNDVI, =8, +Bilnprec;, +6,lnpop, + B,lngdp,, + 

(4) 


考虑 到 面板 数据 可 能 存在 同步 相关 \ 面 板 数据 
异 方差 序列 相关 等 问题 ,在 进行 面板 数据 回归 分 析 
过 程 中 ,选取 面板 矫正 标准 误 (PCSE) 估 计 方 法 对 
定 效 应 模型 进行 修正 ,以 加 强 回归 结果 的 有 效 性 和 
稳健 性 。PCSE 方法 在 保留 最 小 二 乘法 (OLS) 佑 计 
参数 的 基础 上 将 残 差 项 代入 对 角 和 矩阵 ,对 其 标准 差 
进行 修正 ,能 够 更 加 准确 地 对 面板 数据 进行 回归 
估计 。 

变量 解释 及 描述 性 统计 ( 表 1) 如 下 : 


Bufarm， + B;RGPF;, + Mi 


FF 区 
R1 变量 定义 与 基本 统计 
Tab.1 Variable definitions and basic statistics 
变量 变量 描述 均值 标准 差 ”期 望 
InNDVI;。 对 数 变 换 的 NDVI 均值 -0.635 4 0.3690 na 
Inprec,, 对 数 变换 的 年 降水 量 6.388 0 0.3053 + 
对 数 变换 的 人 口 密度 = 人 
1 = 
Npop ix 口 总 数 / 土 地 面积 5.2556 0.739 4 
对 数 变 换 的 经 济 发 展 水 
lngdp， 三 = 
ngdpy。 平 GDP/ 土 地 面积 4.6538 1.0653 
农耕 强度 = 耕地 面积 / 土 
farm, 地 总 面积 0.1893 0.1213 -= 
退耕 强度 = 累计 退耕 还 林 
RGPF, © ; f 
“面积 /土地 总 面积 nee eee 


iE: GDP 经 CPI 调整 (以 1997 年 为 基准 );“ +” (RAS EA A BO 
HE,“ - ”代表 回归 系数 期 望 为 负 。 


(1) 植被 覆盖 水 平 ,用 NDVI 来 表征 。 为 分 析 
退耕 还 林 工 程 对 植被 恢复 的 驱动 效应 ,将 NDVI 设 

(2) 降雨 水 平 ,用 年 降水 量 (prec) 来 表征 。 降 
水 量 对 植物 生长 的 积极 影响 在 相关 研究 中 已 被 充分 
证 实 ,一 般 情 况 下 ,降水 量 越 多 的 地 区 植被 覆盖 水 平 
越 高 ,两 者 成 正 向 关系 。 

(3) 人 地 关系 ,用 单位 面积 人 口 数量 即 人 口 密 
度 (pop) 来 表征 。 其 反映 人 口 数量 与 土地 资源 之 间 
的 关系 。 一 般 情 况 下 ,人口 密 度 越 高 ,土地 需求 和 利 
用 程度 越 高 ,植被 覆盖 水 平 越 低 。 

(4) 经 济 发 展 水 平 。 为 消除 地 区 面积 对 生产 总 
值 的 影响 ,选用 单位 土地 面积 的 GDP ( gdp ) 作为 地 
区 经 济 发 展 水 平 的 表征 指标 。 一 般 情 况 下 ,经 济 发 
展 水 平 越 高 ,对 土地 的 利用 程度 越 高 ,植被 覆盖 水 平 
越 低 。 同 时 ,为 消除 价格 波动 的 影响 ,将 各 市 的 名 义 
国民 生产 总 值 指数 折算 成 1997 年 不 变价 的 实际 人 
均 GDP, 

(5) 农耕 强度 ,用 单位 土地 面积 的 耕地 面积 
(farm) 来 表征 。 其 反映 土地 总 面积 与 耕地 面积 的 数 
量 关系 ,一 般 情况 下 ,农耕 强度 与 植被 覆盖 水 平成 负 
向 关系 。 

(6) 退耕 还 林 工 程 。1999 年 陕西 省 开始 施行 
退耕 还 林 工 程 ,直至 2009 年 第 一 轮 退耕 还 林 工 程 阶 
段 性 结束 。 为 量化 这 一 研究 对 象 , 模 型 选用 单位 土 
地 面积 的 退耕 还 林 面 积 (RGPF ) 作为 研究 指标 。 退 
耕 还 林 工 程 是 在 减少 耕地 面积 的 基础 上 ,人 为 干预 
植物 生长 的 生态 工程 ,因此 ,预期 退耕 还 林 工 程 与 植 
被 覆盖 水 平成 正 向 关系 。 
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3 结果 与 分 析 


3.1 植被 覆盖 年 际 变化 特征 

3.1.1 NDVI 值 年 际 发 展 趋势 ”提取 1998 一 2013 
年 陕西 省 历年 的 NDVI 平 均值 来 表示 NDVI 的 年 际 
变化 特征 (图 2a)。 由 图 2a 可 知 , 研 究 期 内 陕西 省 
NDVI 值 总 体 呈 现 上 升 趋势 ,而 年 际 间 的 NDVI 指数 
存在 波动 ,一 定 程 度 上 是 受降 水 量 等 自然 条 件 的 影 
响 。 其 中 ,2000 年 的 植被 覆盖 在 研究 区 间 内 最 低 ， 
NDVI 年 平均 值 为 0.439 ;2012 年 的 植被 覆盖 在 研究 
期 内 最 高 ,NDVI 年 平均 值 为 0. 601 ;整体 来 看 ,陕西 
省 NDVI 年 平均 值 增长 速率 为 0.95% ,表明 陕西 省 
在 1998 一 2013 年 NDVI 值 不 断 增加 ,植被 恢复 取得 
了 阶段 性 进展 。 

3.1.2 植被 覆盖 度 变化 ”利用 NDVI 数据 分 别 计 
算出 1998—2013 年 陕西 省 10 个 市 的 植被 覆盖 度 
(图 2)。 由 图 2 可 知 ,虽然 各 市 的 植被 覆盖 度 年 际 
间 出 现 波动 ,但 整体 仍 呈 上 升 趋势 。 其 中 ,位 于 陕 北 
退耕 还 林 区 的 延安 市 和 陕 南 退耕 还 林 区 的 安康 市 植 
被 覆盖 度 的 增加 尤为 显著 。 延 安 市 1998 年 的 植被 
覆盖 度 为 39.79% ,至 2013 年 已 增长 为 55.18% ; 安 
康 市 的 植被 覆盖 度 由 1998 年 的 74. 07% 增长 至 
2013 年 的 84.23% 。 值 得 注意 的 是 ,延安 市 地 处 陕 
北 地 区 ,自然 条 件 相 比 其 他 区 域 并 不 占 优 ,但 相 比 其 
他 地 级 市 ,延安 市 植被 覆盖 度 增 长 最 快 ,年 增长 速率 
为 1.12% 。 与 此 同时 ,通过 统计 数据 可 得 ,延安 市 
在 1998 一 2013 年 累计 退耕 还 林 面 积 占 全 省 的 
23.95% ,这 被 视 为 其 植被 覆盖 度 增 加 的 主要 原因 ， 
也 在 一 定 程 度 上 表明 了 退耕 还 林 工 程 对 植被 覆盖 增 
加 存在 积极 影响 。 为 进一步 验证 退耕 还 林 工 程 对 植 
被 覆盖 变化 影响 的 显著 性 ,引入 统计 模型 对 其 进行 
深入 分 析 。 

3.2 植被 覆盖 区 域 差 异 特征 

通过 对 1998 年 .2009 年 与 2013 年 的 SPOT，. 
NDVI 数 据 进行 琶 加 计算 ,可 得 到 年 份 间 的 NDVI 差 
值 图 (图 3)。 由 图 3c 可 知 ,1998 一 2013 年 除 关 中 地 
区 外 ,全 省 大 面积 地 区 植被 覆盖 都 得 到 了 显著 提升 。 
相 较 于 1998 年 ,2013 年 陕西 省 仅 5. 93% 区 域 的 
NDVI 值 轻微 缩减 , 旦 缩减 幅度 小 于 0.2; 而 绝 大 部 
分 区 域 的 NDVI 值 都 得 到 提升 ,其 中 ,40. 61% 区 域 
的 NDVI 值 提升 了 约 0.1,29.05% 区 域 的 NDVI 值 
提升 了 约 0.2,23. 69% 区 域 的 NDVI 指数 提升 了 
0.25 以 上 。 
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图 2 1998—2013 年 陕西 省 NDVI 和 植被 覆盖 度 的 逐年 变化 


Fig.2 Annual changes of NDVI and vegetation coverage in Shaanxi Province from 1998 to 2013 


陕西 省 退耕 还 林 工 程 于 2009 年 开始 进入 验收 
与 维护 阶段 (新 一 期 的 退耕 还 林 工 程 不 在 研究 区 间 
内 , 故 暂 不 研究 ) , 1998—2009 年 及 2009—2013 年 
的 陕西 省 NDVI 值 变化 (图 3a 与 图 3b) 响 应 了 这 一 
人 为 要 素 的 变化 。NDVI 差 值 显示 ,1998 一 2009 年 
陕西 省 大 部 分 区 域 NDVI 指数 都 得 到 较 大 幅度 提 
升 ,而 2009 一 2013 年 只 有 小 部 分 区 域 NDVI 值 发 生 
明显 变化 ,多 数 区 域 基本 维持 稳定 状态 。 经 统计 分 
析 可 知 ,1998 一 2009 年 , 仅 西 安 宝鸡、 渭南 汉中 4 
市 中 的 少数 区 域 NDVI 值 略微 缩减 ,省 内 其 他 区 域 
AY) NDVI 值 均 得 到 不 同 程度 的 增加 ;2009 一 2013 年 ， 
西安 .铜川 宝鸡 3 市 的 NDVI 值 普遍 降低 ,其 他 区 
域 基本 持平 。 这 显然 表明 1998 一 2009 年 植被 状况 
改善 更 为 显著 。 

综 上 所 述 ,1998 一 2013 年 省 内 绝 大 部 分 市 区 
NDVI 值 呈 现 上 升 趋势 ,是 不 同 阶 段 呈 现 出 不 同 的 
变化 特征 。 在 1998 一 2009 年 退耕 还 林 工 程 进行 期 
间 ,研究 区 内 绝 大 部 分 区 域 NDVI 值 呈 一 致 性 增加 
趋势 ,植被 恢复 效果 显著 ,从 侧面 反映 出 退耕 还 林 工 


程 的 实施 对 于 陕西 省 的 植被 恢复 有 着 积极 的 作用 。 
相 比 之 下 ,2009 年 退耕 还 林 工 程 进 入 验收 与 维护 阶 
段 后 ,各 区 域 NDVI 值 变化 增 减 不 一 , 陕 北 地 区 植被 
羡 基 本 不 变 , 陕 南 地 区 植被 覆盖 略微 增加 ,关中 地 

区 各 市 植被 覆盖 普遍 缩减 。 

无 论 从 全 域 还 是 各 市 区 的 植被 覆盖 时 空 变化 来 
看 ,陕西 省 植被 覆盖 水 平均 呈现 改善 趋势 ,但 是 , 仍 
然 不 能 完全 确信 这 一 改善 是 由 退耕 还 林 工 程 实施 推 
动 的 ,也 有 可 能 是 其 他 自然 或 社会 经 济 因素 导致 了 
以 上 观察 结果 。 鉴 于 此 ,笔者 通过 回归 模型 来 定量 
分 析 退 耕 还 林 工 程 实施 对 植被 恢复 的 驱动 效应 。 
3.3 ”退耕 还 林 工 程 对 植被 恢复 的 驱动 效应 

利用 1998 一 2013 年 陕西 省 10 个 市 级 行政 区 的 
16 a 数据 ,构建 面板 数据 模型 来 重点 探讨 退耕 还 林 
工程 实施 与 陕西 省 植被 覆盖 变化 的 相关 关系 ,将 选 
取 的 6 个 因素 值 导 入 市 级 研究 单元 ,实现 因素 的 定 
量化 和 空间 化 ,得 到 其 空间 分 布 图 (图 4)。 

为 了 消除 面板 数据 存在 的 同期 相关 和 截面 异 方 
差 问题 ,选取 PCSE 估 计 方 法 对 固定 效应 模型 进行 
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Fig.3 Difference of values of NDVI in Shaanxi Province in typical years 
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注 :本 文 共 使 用 了 1998—2013 年 16 a 的 数据 ,考虑 到 篇 幅 限 制 , 图 中 仅 选 取 2013 4 
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FE 各 地 市 因素 数据 进行 空间 可 视 化 展示 。 


图 4 陕西 省 植被 恢复 驱动 因素 定量 化 空间 分 布 


Fig.4 Quantitative spatial distribution of driving factors for vegetation regeneration in Shaanxi Province 


修正 ,以 加 强 回归 结果 的 有 效 性 和 稳健 性 ,更 充分 地 
反映 NDVI 与 各 个 可 能 的 影响 因素 之 间 的 关系 。 回 


归结 果 如 表 2 所 示 。 


从 表 2 可 知 , 拟 合 优 度 RR =0.751 5 ,说 明 模型 
整体 拟 合 效果 和 较 好 ,模型 设 定 有 效 。 从 各 变量 对 应 


P 值 可 知 , 在 0.01 检验 水 平 下 ,降水 量 、 人 口 密度 、 
退耕 强度 .地 均 GDP 农耕 强度 变量 对 于 陕西 省 的 


NDVI 均值 有 着 显著 影响 (P <0. 01) , 除 人 口 密度 外 


各 指标 影响 方向 均 与 预期 假设 相符 。 
降水 量 .农耕 强度 .地 均 CDP 人口 密度 4 个 变 


成 


=> 


表 2 PCSE 估计 方法 面板 数据 回归 结果 
Tab.2 Regression results of PCSE 


自 变量 参数 估计 值 标准 差 P 值 
lnprec， 0.459 4 0.091 9 0.00 0*** 
Inpop, 0.513 1 0.030 5 0.000*** 
Ingdp , -0.075 8 0.014 6 0.00 0*** 
farm; -2.743 6 0.189 1 0.000*** 
RGPF,, 1.972 5 0.707 1 0.005 *** 

Bo -0.405 6 0.248 2 0.00 0*** 

调整 后 的 R? 0.751 5 


注 :* 表 示 已 <0.10,** 表 示 已 <0.05,* * * 表 示 已 <0.01。 


量 对 植被 覆盖 水 平 的 影响 各 不 相同 。 从 回归 结果 来 
看 ,降水 量 每 增加 1 个 百分点 ,可 使 NDVI 值 增加 
0.46 个 百分点 。 植 物 生 长 受降 水 量 影响 较 大 ,异常 
的 旱 涝 天 气 都 会 严重 影响 植物 生长 情况 ,生态 脆弱 
的 区 域 对 降水 量变 化 更 敏感 。 农 耕 强 度 每 增加 1 个 
百分点 ,NDVI 值 将 下 降 2.74 个 百分点 ,影响 程度 居 
所 有 变量 之 首 。 粗 放 式 的 农耕 方式 严重 危害 到 陕西 
省 的 生态 环境 ,改善 陕西 省 生态 环境 咪 需 在 满足 粮 
食 产量 需求 的 前 提 下 ,减少 不 适宜 耕作 的 农耕 地 并 
重建 生态 系统 。 经 济 发 展 水 平 每 提升 1 个 百分点 ， 
NDVI 值 将 下 降 0. 075 8 个 百分点 。 经 济 增长 带 来 
的 是 资源 需求 的 增多 和 环境 压力 的 加 大 ,因此 ,在 保 
证 经 济 正常 发 展 的 前 提 下 需 尽 可 能 减 小 其 对 生态 环 
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种 的 耕地 ,栽种 更 适合 生长 的 林 草 ,减少 了 耕地 面 
耻 , 对 植被 恢复 有 着 积极 影响 。 除 上 述 影响 外 ,退耕 
还 林 工 程 还 改变 了 农村 的 生产 .生活 经 营 方式 ,进而 
改善 了 农村 生态 环境 。 随 着 退耕 还 林 工 程 的 逐步 推 
进 ,农村 大 量 劳 动力 从 粮食 生产 中 解放 出 来 ,从 事 种 
植 业 养殖 业 ,设施 农业 等 劳务 经 济 活动 ,改变 了 农 
村 的 现 有 产业 结构 ; 随 之 而 来 的 是 退耕 农户 收入 的 
增长 以 及 生活 方式 的 改变 ,使 得 农村 户 用 沼气 得 到 
迅速 发 展 , 煤 太阳能、 天 然 气 等 逐渐 取代 烧 柴 ,环境 
压力 得 到 释放 。 从 长 远 角 度 来 看 ,退耕 还 林 工 程 在 
改善 生态 环境 的 同时 ,更 增强 了 全 社会 的 生态 保护 
意识 ,是 生态 环境 持续 改善 的 根本 保障 。 


4 结论 


(1) 从 NDVI 时 序 变 化 结果 来 看 ,1998 一 2013 
年 陕西 省 植被 覆盖 水 平 总 体 呈 上 升 趋势 ,NDVI 值 
年 平均 增长 速率 为 0.95% ,植被 恢复 效果 明显 。 通 
过 对 比特 定年 份 间 的 NDVI 差 值 ,可 以 看 出 , 相 比 
2009 一 2013 年 ,1998 一 2009 年 ( 即 陕西 省 全 面 实施 
退耕 还 林 工 程 期 间 ) 研 究 区 NDVI 值 提 升 尤为 显著 。 

(2) 从 植被 覆盖 度 空间 分 布 情况 来 看 ,陕西 省 
植被 覆盖 水 平整 体 上 由 南 向 北 递 减 . 空 间 分 异 特征 


境 的 影响 ,实现 可 持续 发 展 。 与 预期 不 同 的 是 ,人 口 
密度 的 增加 对 NDVI 值 有 着 正面 的 影响 。 这 一 方向 
与 预期 的 偏差 在 一 定 程 度 上 是 陕西 省 城镇 化 率 的 快 
速 提升 导致 的 。 近 年 来 ,陕西 省 城镇 化 步伐 不 断 


显著 。16 a 间 ,省 域内 NDVI 值 有 不 同 程度 的 上 升 ， 
40. 61% 的 区 域 提升 约 0.1,29. 05% 的 区 域 提 升 约 
0.2,23.69% 的 区 域 提 升 0.25 以 上 。 分 析 表 明 在 市 
域 尺度 上 植被 覆盖 的 增加 与 退耕 还 林 工 程 相关 性 明 
显 ,植被 恢复 效果 最 好 的 两 个 市 (区 ) 分 别 为 延安 市 


加 快 ,1998 一 2010 年 ,城镇 人 口 比 重 整 体 呈 增长 趋 
势 ,由 21.80% 增加 到 50.02% 。 城 镇 化 程度 越 高 ， 
人 口 越 集中 , 非 城镇 区 域 的 环境 压力 随 之 降低 。 
值得 注意 的 是 ,退耕 还 林 工 程 的 实施 对 陕西 省 
植被 恢复 具有 明显 正面 效应 。 由 表 2 可 以 看 出 , 退 
耕 还 林 面 积 占 总 面积 的 比值 每 提升 1 个 百分点 , 预 
计 可 使 NDVI 值 增加 1. 97 个 百分点 ,换算 成 降水 量 
来 看 ,提升 陕西 省 退耕 还 林 面 积 占 总 面积 比值 1 个 
百分点 ,相当 于 陕西 省 的 降水 量 增加 了 4. 29 个 百 分 
点 。 这 一 结果 对 降水 量 偏 少 的 西部 地 区 具有 重要 意 
义 。 退 耕 还 林 工 程 的 实施 可 增加 植被 覆盖 度 ,提升 
生态 系统 的 抵御 能 力 , 减 小 旱 涝 等 气象 灾害 对 植被 
造成 的 负面 影响 。1999 一 2013 年 ,陕西 省 累计 完成 
退耕 还 林 2.41 x 10 hm ,其 中 退耕 地 还 林 1. 02 x 
10° hm? ,荒山 造林 1. 26 x 10° hm? ,封山育林 1. 29 x 
10° hm ”。 退 耕 还 林 工 程 退 出 坡地 及 其 他 不 适宜 耕 


〈 陕 北 退耕 还 林 区 ) 、 安 康 市 ( 陕 南 退耕 还 林 区 ) ,其 
中 延安 市 的 植被 覆盖 水 平 增幅 最 大 ,年 增长 率 达 
1.12%, 

(3) 退耕 还 林 工 程 的 实施 对 陕西 省 植被 恢复 具 
有 明显 的 积极 效应 。 通 过 统计 陕西 省 1999 一 2009 
年 累计 退耕 还 林 面 积 ,并 引入 降雨 水 平 人 地 关系 、 
经 济 发 展 水 平 .农耕 强度 等 自然 .社会 经 济 因 素 变 
量 , 与 对 应 地 区 8 月 的 NDVI 平均 值 进行 回归 分 析 。 
结果 显示 ,退耕 还 林 面 积 占 总 面积 的 比值 每 提升 1 
个 百分点 ,预计 可 使 NDVI 值 增加 1. 97 个 百分点 。 
这 一 定量 分 析 进 一 步 验 证 了 时 间 动 态 分 析 的 结 
说 明 退 耕 还 林 工 程 对 陕西 省 植被 恢复 的 影响 是 积极 
且 显 著 的 。 

需要 指出 的 是 ,植物 生长 是 多 种 因素 综合 作用 
的 结果 , 除 所 列举 的 因素 外 ,不 同 区域 的 土壤 状况 、 


FOR 区 


地 形 地 貌 政策 条 件 .产业 结构 等 因素 对 该 区 域 植被 

盖 变 化 同样 起 着 关键 作用 。 以 榆林 、 延 安 两 市 为 
例 , 二 者 同属 陕 北 地 区 ,但 榆林 地 理 位 置 更 为 偏 北 ， 
且 二 者 产业 结构 不 尽 相同 ,虽然 榆林 累计 退耕 还 林 
面积 在 全 省 范围 内 仪 次 于 延安 ,但 其 植被 恢复 效果 
相 比 延安 市 差距 较 大 。 鉴 于 此 ,后 续 研 究 应 考虑 加 
入 地 形 地 貌 .产业 结构 等 相关 因素 。 与 此 同时 ,本文 
使 用 了 遥感 数据 分 析 植 被 恢复 状况 , 受 遥 感 数据 分 
辩 率 限制 , 解 译 精度 受到 一 定 影响 。 在 今后 的 研究 
中 , 若 采 用 更 高 空间 分 辩 率 的 遥感 影像 ,将 可 以 更 加 
准确 客观 地 评价 植被 恢复 动态 ,进一步 深入 分 析 植 
被 恢复 过 程 中 出 现 的 问题 。 值 得 注意 的 是 ,陕西 省 
地 域 差异 大 ,各 地 自然 环境 、 经 济 结构 不 尽 相 同 ,时 
空 特征 所 展现 出 的 是 全 省 及 省 内 不 同市 区 植被 覆盖 
情况 对 退耕 还 林 工 程 的 响应 。 这 表明 ,从 大 尺度 研 
究 单元 上 可 区 分 出 植被 覆盖 变化 过 程 中 生态 工程 的 
信号 。 但 随 着 生态 文明 建设 的 深入 推进 ,对 具体 微 
观 尺 度 (如 样 带 、 县 域 镇 域 等 ) 植 被 恢复 建设 工作 
进行 科学 合理 的 指导 显得 尤为 迫切 ,因此 ,今后 仍 需 
要 在 微观 尺度 上 对 相关 问题 进行 深入 探讨 。 
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Effect of Grain for Green Project on Vegetation Regeneration in 
Shaanxi Province : Based on PCSE Panel Data Model 
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Abstract; In this paper, the temporal and spatial variations of vegetation cover in Shaanxi Province were ana- 
lyzed. The purposes of this study were to explore the driving effect of the Grain for Green Project on vegetation re- 
generation, and the relative contributions of some natural and socio-economic factors (such as the Grain for Green 
Project and rainfall level) to vegetation cover change were estimated based on the SPOT), NDVI data during the 
period from 1998 to 2013 by using the pixel dichotomy method and the Panel Corrected Standard Errors (PCSE ) . 
The result showed that the vegetation coverage in Shaanxi Province was holistically improved, the annual average 
NDVI value was increased by 0.98% from 1998 to 2013, and the effect of vegetation regeneration was obvious by 
implementing the Grain for Green Project from 1998 to 2009. The spatial difference of vegetation regeneration effect 
in the province was significant. During the 16 years,the NDVI in the province increased with different rates, it was 
increased by 0.1 in 40.61% area,by 0.2 in 29.05% area,and by 0.25 in 23.69% area. In which the most evi- 
dent effect of the project on vegetation regeneration occurred in Yan’an City located in north Shaanxi Province , 
where the NDVI was increased by 1. 12% . The NDVI was increased by 1.97% if the Grain for Green Project was 
implemented in 1% area of the total, which was equivalent to increasing 4.29% precipitation. The results of this 
study are of great significance for reconstructing the ecological environment in the study area. 
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